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■Abstragt 

Axisymmetric subsonic diffuser flows were calcu- 
lated with the NPARC Navier-Stokes code in order to 
determine the effects various code features have on 
the flow solutions. The code features examined in 
this work were turbulence models and boundary con- 
ditions. Four turbulence models available in NPARC 
were used: the Baldwin-Lomax algebraic model, the 
Baldwin-Barth one-equation model, and the Chien k- 
e and Wilcox k -(0 two-equation models. The three 
boundary conditions examined were the free bound- 
ary, the mass flux boundary and the subsonic outflow 
with variable static pressure. In addition to boundary 
condition type, the geometry downstream of the dif- 
fuser was varied to see if upstream influences were 
present. The NPARC results are compared with 
experimental data and recommendations are given for 
using NPARC to compute similar flows. 

Introduction 

The NPARC Navier-Stokes code^ is used by gov- 
ernment, industry and academia to calculate a wide 
variety of aerospace propulsion flows. NPARC is 
currently being used to calculate subsonic diffuser 
flows encountered in both subsonic and supersonic 
aircraft engine inlets. These flows are turbulent and 
characterized by strong adverse pressure gradients. 

Predicting turbulent adverse pressure gradient 
flows, both with and without separation, is a chal- 
lenging task for most turbulence models and flow 
solvers. Previous work has shown, for example, that 
the Baldwin-Barth one-equation model^ and the 
Baldwin-Lomax algebraic model^ generally under 
predict the shear stress for these flows, while k-£ 
models tend to over predict wall shear stress. In 
flows with separation or incipient separation, this 
results in the Baldwin-Barth and the Baldwin-Lomax 
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models predicting early separation, and k-e predicting 
late separation. 

Subsonic diffuser studies with the NPARC code, 
described herein, reinforce these turbulence model 
observations and emphasize the need for the capabil- 
ity to calculate adverse pressure gradient flows. In 
response to this, the k-co model of Wilcox,^ which is 
known to give better results for adverse pressure gra- 
dient flows, has recently been installed in NPARC.® 

The objective of this paper is to evaluate NPARC 
for calculating subsonic diffuser flows, with emphasis 
on the effects that turbulence model and boundary 
condition selection have on the quality of the flow 
solution. Two different experimental diffuser geome- 
tries are examined. The first is a diffusing pipe flow 
referred to as Fraser Flow A in the AFOSR-IFP-Stan- 
ford Conference Proceedings;^’^^ this flow remains 
attached. The second geometry consists of a turbu- 
lent boundary layer developing axially on a cylinder 
in a wind tunnel with diverging walls. Two dif- 
ferent pressure gradient flows were measured for this 
case: one which remains attached and one which is 
separated. 

The following sections provide a brief description 
of the turbulence models in NPARC and a description 
of how NPARC was used to calculate each of the dif- 
fuser flows. The computational results are evaluated 
based on their agreement with experimental data. 
The flows examined were axisymmetric and were run 
using the two-dimensional NPARC code. NPARC 
version 2.1 was used for all calculations except those 
using the k-(D turbulence model; these cases were run 
using NPARC version 2.2. 


Tlirbulence Models in NPARC 
The turbulence models used in this study were the 
Baldwin-Lomax algebraic model (BL),® the Baldwin- 
Barth one-equation model (BB),^ the Chien low Rey- 
nolds number k-e model^^’*'^ and the recently 
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installed Wilcox k-co model^’* (both two-equation mod- 
els). 

The Baldwin-Lomax algebraic model is patterned after 
the model of Cebeci'^ with modifications that avoid the 
necessity for finding the boundary layer edge. It is a 
two-layer model which makes use of the Prandtl-Van 
Driest formulation for the inner layer.^ 

The Baldwin-Barth one-equation model avoids the 
need for an algebraic length scale and is derived from a 
simplified form of the k-e equations. It solves a field 
equation for the turbulence Reynolds number, Rx=kVve. 

The Chien low Reynolds number k-e model solves two 
transport equations - one for the turbulent kinetic 
energy, k, and one for the turbulent dissipation rate, e - 
with the turbulent viscosity proportional to k^/e. The 
phrase, “low Reynolds number,” refers to the fact that 
the model is applied near the solid surface where the tur- 
bulent Reynolds number and wall functions are not 
required. 

The Wilcox k-co model also solves two transport equa- 
tions, where the second quantity, co, is the dissipation 
divided by the turbulent kinetic energy. This model has 
been shown to predict adverse pressure gradient flows, 
both with and without separation, better than the other 
three turbulence models,^’ 

Fraser Diffusing Pipe Flow 

The first flow case that was calculated is known as the 
Fraser (flow A) case from the AFOSR-IFP Stanford Con- 
ference.^’^^ The geometry of the conical diffuser is 
shown in figure 1. In the experiment, a length of 0.152 
m diameter straight pipe preceded the 5 degree half 
angle conical diffuser. The core flow velocity at the dif- 
fuser entrance was approximately 52 m/s (Mach 0.15) 
and Reg, the Reynolds number based on the momentum 
thickness 0, was approximately 3000. Measurements of 
velocity profiles and skin friction were made at 11 loca- 
tions in the diffuser as listed in table 1. The first mea- 
surement station, corresponding to x=0.117 m (diameter 
= 0.154 m), is slightly upstream of the beginning of the 
diffusing section; the last measurement station, corre- 
sponding to x=0.642 m (diameter = 0.236 m), is just 
upstream of the exit plane of the diffusing section. The 
diffuser exits into ambient air. 
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Figure 1 . Schematic of Fraser Flow A experimental setup. 


Table 1 . Locations of velocity measurements 
for Fraser Flow A. 


Station Number 

Axial Location (m) 

1 

0.117 

2 

0.158 

3 

0.211 

4 

0.236 

5 

0.300 

6 

0.381 

7 

0.438 

8 

0.468 

9 

0.522 

10 

0.578 

11 

0.642 


Computational Grids 

The baseline computational grid for this case was gen- 
erated using a version of the INGRID grid generation 
package*^ and has 121 points in the axial direction and 
71 points in the radial direction. As shown in figure 2, 
the grid includes a straight section of pipe both upstream 
and downstream of the conical diffusing section. A short 
circular arc transition was used between sections. The 
grid was packed to the solid wall such that the value of 
y+ at the first point off the wall was approximately 1; 
this value is based on the core flow velocity at the first 
measurement station and an assumed skin friction coeffi- 
cient of 0.003. Previous NPARC validation work indi- 
cates that this gives sufficient boundary layer resolution 
for wall bounded turbulent flows. The grid is also 
packed axially at the inflow boundary, as recommended 
in reference 17, to resolve the large axial flow gradients. 



Figure 2. Baseline computational grid for Fraser Flow A. 


Grids F2 and F3, shown in figure 3, were also gener- 
ated, in addition to the baseline grid. These grids were 
different from the baseline grid primarily in the geome- 
try downstream of the diffusing section; they were of 
interest because it was initially suspected that the down- 
stream geometry had some upstream influence on the 
flow solutions. Grid F2 is identical to the baseline grid, 
except that the dimensions are 181x71, and the axial grid 
points are clustered both at the inflow boundary and at 
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the last measurement station in the diffuser. Grid F3 has 
the same dimensions as grid F2, but has different pack- 
ing and a different geometry downstream of the last mea- 
surement station. It is packed at the inflow boundary, 
the last measurement station and the outflow boundary. 
Instead of using a straight section downstreeun of the last 
measurement station, grid F3 uses an extension of the 5 
degree half angle conical diffusing section. 



(a) Grid F2 



(b) Grid F3 

Figure 3. Alternate computational grids for Fraser Flow A. 


Boundary Conditions 

For the boundary conditions, the inflow of the calcula- 
tions was specified as a free boundary and placed 0.61 m 
upstream of the first measurement station in an attempt 
to match the momentum thickness and displacement 
thickness measured in the experiment at this first station. 
A free inflow boundary in NPARC requires the total 
pressure and temperature to be specified, and uses simple 
characteristic equations to specify the flow field. The 
pipe centerline boundary was specified as an axis of 
symmetry, and the solid boundary was specified as a no- 
slip adiabatic surface. For calculations made using the 
baseline grid, the no-slip, adiabatic boundary condition 
was set for the entire axial length of the pipe. For calcu- 
lations made using grids F2 and F3, the no-slip, adiabatic 
condition was set from the inflow boundary to the last 
measurement station in the diffuser, and a slip surface 
was set for the remaining axial length. 

Using the Chien k-£ turbulence model, the outflow 
location, geometry and boundary condition type were 
examined to determine their effect on the flow solution. 
Using the baseline grid and the specified mass flux 
boundary condition, three outflow locations were exam- 
ined: (1) the last measurement station (the true diffuser 


exit plane), (2) 0.15 m downstream of the last station, 
and (3) 0.37 m downstream of the last station. These 
three locations are labeled in figure 2. The latter two 
outflows were examined because the flow solution 
obtained with the exit at the last measurement station did 
not conserve mass between the outflow grid line and the 
grid line just upstream. The two cases with the extended 
grid conserved mass and produced essentially the same 
solution for the quantities examined: the local skin fric- 
tion coefficient, the static pressure coefficient, the inte- 
gral boundary layer properties, and the velocity 
profiles. With the outflow at location 2, the solutions 
obtained using the baseline grid and grids F2 and F3 
were compared and found to have no appreciable differ- 
ences. 

Using the baseline grid, three different types of exit 
boundary conditions were examined: the free boundary, 
the subsonic outflow boundary with variable static pres- 
sure, and the specified mass flux boundary. These three 
boundary conditions all are obtained using extrapolation 
of the upstream flow field, but use different static pres- 
sure profiles. The free boundary uses the user-specified 
value of the static pressure across the entire boundary; 
the subsonic outflow with variable static pressure uses 
the user-specified value of static pressure at a specified 
grid point with the pressure variation across the bound- 
ary coming from the closest upstream station; and the 
specified mass flux boundary uses the user-specified 
mass flux to compute the static pressure which is then 
held constant across the entire boundary. All three 
boundary conditions produced essentially the same flow 
solution. For the results which follow, the exit boundary 
is located at x=0.15 m downstream of the last station and 
the specified mass flux boundary condition was used. 

This parametric investigation indicates that the flow 
solution for the Fraser subsonic diffuser is essentially 
independent of the outflow geometry, location and 
boundary condition type. The results which follow were 
computed using the baseline grid with the downstream 
boundary at location 2, and using the specified mass flux 
outflow boundary condition. 

Artificial Viscositv 

The artificial viscosity selected in NPARC for these 
computations was the modified Jameson-style artificial 
viscosity with the Jameson-style spectral radius term. To 
remove any doubts that the artificial viscosity was influ- 
encing the turbulent viscous effects, the second-order 
artificial viscosity coefficient was set to zero, although 
increasing it to 0.25 (the default value) had very little 
effect on the flow solutions. The fourth-order coefficient 
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was left at its default value of 0.64, as recommended in 
the NPARC User’s Guide.* 

Comparison of NPARC Solutions to Experimental Data 
The computational results are compared with the 
experimental piezometer and pitot measurements in fig- 
ures 4 through 7. The integral momentum and displace- 
ment thicknesses, 0 and 6*, shown in figure 4a and 4b, 
agree well with the experimental values at the first mea- 
surement station, indicating that the upstream boundary 
conditions were adequately specified. Further down- 
stream in the diffuser, agreement is fairly good with the 
k-e, k-o>, and BL models only slightly under predicting 
the displacement effects, and the BB model over predict- 
ing them. A plot of the shape factor, H, is shown in fig- 
ure 4c, and indicates that all four turbulence models have 
difficulty predicting the adverse pressure gradient 
effects, with the BB and BL models predicting stronger 
retardation effects, while both two-equation models 
under predict these effects. This behavior is consistent 
with the skin friction results of figure 5, which indicate 
that BB predicts separation at x=0.323 m and BL pre- 
dicts separation at x=0,438 m. The k-e skin friction is 
much too high, whereas the k-m results are in very close 
agreement with the experimental data. The pressure dis- 
tribution is shown in figure 6, where 

== ^^P“Pref^'^vPrcf‘^rcfJ ’ reference pressure, 

density and velocity, and are the values at 

station 1, just upstream of the diffusing section. In the 
downstream portion of the duct, BB under predicts the 
pressure rise while the other three turbulence models 
slightly over predict it; this behavior is consistent with 
the displacement effects shown in figure 4b. Three 
velocity profiles are shown in figure 7 corresponding to 
experimental stations 1, 6 and 11 at x = 0.117 m, 0.381 m 
and 0.642 m. All of the models show good agreement at 
station 1, with the exception of the BB model, which 
overshoots the experimental data near the edge of the 
boundary layer. Similar overshoots are also present in 
the NPARC validation of turbulent flow over a flat 
plate, and to a much lesser extent in the adverse pres- 
sure gradient computations of Menter."* 



X (meters) 

(a) Momentum thickness. 



X (meters) 

(b) Displacement thickness. 



X (meters) 


(c) Shape factor. 

Figure 4. Integral boundary layer parameters for Fraser 
Flow A. 
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Figure 5. Local skin friction coefficient for Fraser flow A, 
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(b)x = 0.381m 
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Figure 6. Static pressure coefficient for Fraser flow A, 



0.0 0.2 0.4 0.6 0.8 1.0 1.2 
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(a) x = 0.1 17 m 



U/Ue, 

(c) X = 0.642 m 

Figure 7. Velocity profiles for Fraser flow A at three axial 
locations. 

Driver Adverse Pressure Gradient Flow 
This validation case models the experiments of 
Driver,^^’ who measured flow over a 0.140 m diameter 
cylinder, mounted axially in a wind tunnel with diverg- 
ing walls. A schematic of the experimental setup is 
shown in figure 8. Two types of flow conditions, differ- 
ing in the strength of the imposed pressure gradient, 
were examined. In case BSO, the flow remains attached, 
and in case CSO, the pressure gradient is severe enough 
to cause separation. The pressure gradient was imposed 
on the downstream portion of the flow by diverging the 
four tunnel walls; boundary layer suction was used to 
prevent separation on these walls. Detailed three-dimen- 
sional measurements were used to validate the axisym- 
metry of the flow. Boundary layer profiles were 
measured at 10 stations for case BSO, and 13 stations for 
case CSO, as shown in table 2, where the imposed 


5 





adverse pressure gradient region begins at an axial loca- 
tion of approximately -0.304 m, and the Reynolds num- 
ber based on momentum thickness, Re 0 , was 
approximately 3000. The tunnel was operated at atmo- 
spheric temperature, and the core flow velocity at the 
reference station, station 1 at x=-0.457 m, was 30 m/s 
(M=0.087). 



Figure 8 . Schematic of Driver experimental setup. 


Table 2. Locations of skin friction and velocity 
measurements for the Driver adverse 
pressure gradient experiments. 


Case BSO 

Case CSO 

Station 

Number 

Axial 

Location 

(m) 

Station 

Number 

Axial 

Location 

(m) 

1 

-0.457 

1 

-0.457 

2 

-0.330 

2 

-0.330 

3 

-0.229 

3 

-0.229 

4 

-0.152 

4 

-0.152 

5 

-0.076 

5 

-0.076 

6 

-0.013 

6 

-0.013 

7 

0.013 

7 

0.013 

8 

0.152 

8 

0.025 

9 

0.229 

9 

0.051 

10 

0.305 

10 

0.102 



11 

0.152 



12 

0.229 



13 

0.305 


Computational Grids 

The computational grids for this case, one for case 
BSO and one for case CSO, were also generated using the 
INGRID grid generation package and had 181 points in 
the axial direction and 71 points in the radial direction. 
These grids also included a straight section of pipe both 
upstream and downstream of the diffusing sections, as 
shown in figure 9. The upstream length of pipe, 0.6 m in 


length, was chosen in order to match the momentum 
thickness at the first measurement station; the down- 
stream length of pipe, 0.47 m in length, was chosen to be 
similar in proportion to the downstream length of pipe in 
the Fraser baseline grid. To define the inviscid bound- 
ary, opposite the cylinder surface, the inviscid stream- 
lines, available from the experimental data, were used. 
The intended use of this experimental streamline data 
was for definition of the far field boundary for axisym- 
metric computations. This was needed because the flow 
field around the cylinder was axisymmetric, though the 
test section was not. These streamlines were a distance 
of approximately one to two boundary layer thicknesses 
from the cylinder, so they were translated to increase this 
distance to approximately four to five boundary layer 
thickness, which was more desirable in NPARC. The 
velocity measurements at each station were used when 
translating the streamlines to insure that the same 
amount of mass flux was added at each station. These 
new translated streamlines were then curve fitted, using 
cubic spline interpolation, to give a smooth curve. 
Streamline data was available at only 10 stations for case 
BSO and 13 stations for case CSO. There is some ques- 
tion as to whether or not this amount of data was suffi- 
cient to accurately define the test geometry in the 
computational grid. The grid was packed both axially at 
the inflow boundary and radially at the cylinder surface, 
as described for the Fraser baseline grid. 



(a) Case BSO 



(b) Case CSO 

Figure 9. Computational grids for Driver cases BSO and 
CSO. 
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Boundary Conditions 

The boundary conditions for this case were similar to 
those set for the Fraser case, except for the inviscid 
streamline boundary, which was set to be a slip surface. 

Artificial Viscosity 

The second order artificial yiscosity coefficient was set 
to 0.01, and the fourth-order coefficient, was 0.64. 

Comparison of NPARC Solutions to Experimental Data 

The computational results are compared with experi- 
mental laser-doppler yelocimeter and oil-flow laser 
interferometer data in figures 10 through 13 for case BSO 
and figures 14 through 17 for case CSO. The quantitatiye 
relationships between turbulence models are similar to 
the Fraser flow A results; however, the agreement with 
the experimental data differs somewhat. 

For case BSO, the integral momentum and displace- 
ment thicknesses are shown in figure 10; the agreement 
with experimental data is good in the upstream portion 
of the diffuser. However, both quantities are under pre- 
dicted in the downstream portion of the diffuser, indicat- 
ing an under prediction of the strength of the adverse 
pressure gradient. BB did a better job of predicting the 
adverse pressure gradient effects, as shown by the shape 
factor, figure 10, than the other three models which all 
showed significantly lower H values. The skin friction 
results, shown in figure 11, indicate that all models pre- 
dict attached flow in the duct, with k-8 predicting the 
highest shear stresses, followed by k-co, BL, and BB, 
which is in closest agreement to the data. All of the 
models under predict Cp in the adverse pressure gradient 
region of the duct by approximately the same amount, as 
shown in figure 12. Three velocity profiles are shown in 
figure 13 at experimental stations 1, 5 and 10, at x=- 
0.457 m, -0.076 m, and 0.305 m. All models give fairly 
good agreement with the data at the first experimental 
station, with the agreement degrading further down- 
stream in the duct. 



X (meters) 

(a) Momentum thickness. 



X (meters) 

(b) Displacement thickness. 
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(c) Shape factor. 

Figure 10. Integral boundary layer parameters for Driver 
case BSO. 
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Distance From Wall, meters 



X (meters) 

Figure 11. Lx)cal skin friction coefficient for Driver case 
BSO. 



X (meters) 

Figure 12. Static pressure coefficient for Driver case BSO. 
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(a) X = -0.457 m 



0.0 0.2 0.4 0.6 0.8 1.0 1.2 


u/u^ 

(b) X = -0.076 m 



Figure 13. Velocity profiles for Driver case BSO at three 
axial stations. 

In case CSO, the displacement and momentum thick- 
nesses, shown in figure 14, were also under predicted in 
the downstream portion of the duct, with BB predicting 
the highest displacement thickness and therefore the 
largest response to the adverse pressure gradient. All 
models do a better job of predicting the momentum 
thickness than they did for case BSO. Figure 14c shows 
that all models under predict the shape factor with BB 
giving the highest values, indicating separated flow; the 
other models all give fairly flat profiles. This is consis- 
tent with the skin friction data shown in figure 15, in 
which BB is the only model which predicts separation, 
with a separation bubble beginning at x=-0.079 m and 
ending at x=s0.286 m. In the experiment, the flow 
detached in the vicinity of x=-0.030 m and reattached in 
the vicinity of x= 0.200 m. The trend in the local skin 
friction coefficient is the same as the other two cases. 
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with BB predicting the lowest skin friction, and BL, k-co 
and k-£ predicting successively higher values. The pres- 
sure coefficient, shown in figure 16, was predicted fairly 
well by all of the turbulence models, and significantly 
better than the case BSO results. The velocity profiles at 
stations 1, 5 and 13 are shown in figure 17, and the 
agreement with experimental data is similar to the results 
for cases BSO and Fraser flow A. 



X (meters) 

(a) Momentum thickness. 



X (meters) 

(b) Displacement thickness. 



X (meters) 

(c) Shape factor. 

Figure 14. Integral boundary layer parameters for Driver 
case CSO. 



X (meters) 

Figure 15. Local skin friction coefficient for Driver case 
CSO. 



X (meters) 

Figure 16. Static pressure coefficient for Driver case CSO. 
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(a) X = -0.457 m 
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Figure 17. Velocity profiles for Driver case CSO at three 
axial stations. 


Di.scussion of Results 

For the Fraser flow A conical diffuser case, three dif- 
ferent outflow boundary conditions were used: the free 
boundary, the subsonic outflow boundary with variable 
static pressure, and the specified mass flux boundary. 
All gave valid answers, and may be considered appropri- 
ate for use with subsonic diffuser flow. The geometry 
downstream of the diffusing section was varied as well. 
All geometries worked well, with the exception of the 
baseline grid ending at the diffuser exit; this grid did not 
conserve mass near the exit boundary. 

The Fraser and Driver test cases were run using the k- 
£, k-co, BB, and BL models. The Fraser solutions were 
similar to the results of Menter,"^ but the Driver results 
showed some differences. The trends were similar for 
the displacement thickness, skin friction, and velocity. 
The results differed for the momentum thickness, shape 
factor and static pressure coefficient. The skin friction 
results indicate that NPARC’s ability to predict separa- 
tion is fairly unreliable. For example, the k-£ model pre- 
dicted shear stresses which are too high, and therefore it 
will be likely to predict attached flow when the flow is 
actually separated. The k-to model also tended to predict 
high values of the skin friction coefficient; however, in 
the Fraser case, it was in best agreement with the experi- 
mental data, whereas in the Driver case, it failed to pre- 
dict the separation for case CSO. The BB and BL models 
bracketed the experimental skin friction values, both 
incorrectly predicting separation for the Fraser flow, and 
predicting high values for the Driver BSO case. Of all of 
the models, BB was in the best agreement with the 
Driver experimental data; it was the only model which 
accurately predicted separation for the Driver CSO case, 
while all the other models predicted overly high values 
of shear stress. 

In evaluating these results, some of the factors leading 
to uncertainties should be mentioned. Firstly, the differ- 
ences in experimental measurement techniques should 
be noted. The Fraser data was taken with well-proven 
but less sophisticated measuring equipment (pitot tubes 
and piezometers) than was used to take the Driver data 
(LDV). Another consideration is the quality of the com- 
putational simulation of the experimental conditions. 
The diffusing section for the Fraser flow case was well 
defined, a straight 5 degree half-angle conical diffuser^ 
whereas the Driver diffuser geometry had more uncer- 
tainty because it was defined from the experimental 
streamlines available only at the defined measurement 
stations. Ambiguities in the use of this data, for example 
in the quality and smoothness of the curve fit, may have 
led to the simulation of a pressure gradient different 
from the experimentally imposed gradient. These fac- 
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tors bring to light the possibility that the inconsistencies 
in the validation results may be partially due to the qual- 
ity and usage of the experimental data. 

Conclusions 

The NPARC Navier-Stokes code was applied to axi- 
symmetric subsonic diffuser flows in order to evaluate 
the effects of various code features. The emphasis of 
the study was on outflow boundary condition type, the 
geometry downstream of the diffusing section, and the 
turbulence model used. The results showed that all three 
outflow boundary conditions - the free boundary, the 
specified mass flux boundary, and the subsonic outflow 
boundary with variable static pressure - work well for 
this type of flow problem. All of the downstream geom- 
etries examined also worked equally well, with the 
exception of the Fraser flow baseline grid which ends at 
the last measurement station. The four turbulence mod- 
els available in NPARC - the Baldwin-Lomax (algebraic 
model), the Baldwin-Barth (one equation model), and 
the Chien k-e and Wilcox k-OD (two equation models) - 
were all examined and gave less conclusive results, with 
the k-e model over-predicting the shear stresses in all 
cases, and the other models varying in performance. 
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